THE MOMENT GENERATING FUNCTION FOR RAY LENGTHS IN THE 
HALF GILBERT MODEL WITH RECTANGULAR CELLS 
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| Abstract. In the full rectangular version of Gilbert's tessellation [HO [3] lines extend either 

horizontally (with east- and west-growing rays) or vertically (north- and south-growing rays) 
from seed points which form a Poisson point process, each ray stopping when another ray is met. 
In the half rectangular version [3J, east and south growing rays do not interact with west and 
north rays. Using techniques developed in our previous paper [3j, we derive an exact expression 
^ C"| for the moment generating function for the ray length distribution in the half rectangular model. 
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1. A BRIEF REVIEW OF STOPPING SETS IN THE HALF GILBERT MODEL 

Suppose that a stationary Poisson process of intensity A exists in the plane, with seeds marked 
either H (east growing) with probability q or V (south growing) with probability 1 — q. These 
seeds produce rays growing at a constant rate in directions matching their label (H or V). Rays 
end their growth at the instant that the growing tip meets another ray. Our aim is to investigate 
| the terminal length distribution of an east growing test ray, whose seed location we take to be 

the origin O. The seeds that are relevant for the test ray lie in the unbounded octant between 
If) ', the lines y = x and y = 0, with x > 0: we call this region, the initial live zone. 

We construct a stopping set process [H [5] by expanding a domain — an isosceles right angle 
triangle (see Figure [I]) — into the live zone, stopping when it hits the first seed s\ whose 
coordinates relative to O are (x\,yi). This creates a domain S\ with area E\ that is exponentially 
distributed. If s\ is V-type, then it will provide the ray that blocks the test seed; thus L = x\ 
^ . and no other seeds need be considered. 

Alternatively if s\ is H— type, then, instead of growing S\ (retaining its shape as an isosceles 
right-angle triangle), we remove a part of the live zone: a 'dead zone' labelled D\ (see Figure [1]) 
which has now become irrelevant. As S\ U D\ has been constructed without drawing upon any 
information taken from outside S\ UDi, the point process in the remaining region (the new live 
zone) is still a Poisson process with unchanged intensity given the information within S\ U D\. 

We now grow a trapezium whose left-hand side located at x = x\ has length y = y\. The 
trapezium expands until its right-hand side first hits a seed S2 (in the new live zone). The 
stopping set formed is called S^- It has an exponentially distributed area Ei- 

We proceed in this way, forming a sequence of stopping sets (illustrated in Figure [T]) with in- 
dependent exponentially-distributed areas. The sequence possesses a recursive structure, which 
we can exploit. It is also important to note that the first V-type seed will provide the ray which 
blocks the test ray. Without our introduction of dead zones, a complicated algorithm would be 
required to check if a y-type ray actually reaches the path of the test ray. 

The recursive structure commencing with a generic live zone: Suppose that we begin 
observing the process when the live zone has left boundary of height y and when we are about 
to construct S n . In Figure [H we draw the case n = 2. The probability density function for 



Date: May 7, 2013. 



1 



2 



JAMES BURRIDGE*, RICHARD COWAN** 




the length, r, of S n 7 s base, conditional on the height y of its left boundary, follows from the 
exponential distribution of 5 n 's area E n . It is therefore: 

f(r | y) = \(r + y)e-^ r2+2ry \ 

If the stopping seed s n for set S n is V-type, then its south ray will be the first to intersect the 
test ray and the process ends. Otherwise, another dead zone is created and further trapezoidal 
stopping sets are formed until a V-type is met. 

Let X be the random variable equal to the horizontal distance covered by stopping sets until 
the process comes to an end. The density function of X, conditional on y will be: 



foe r pr-\-y 

g(x\y) = (1 -q)\{x + y)e-^ x2+2xy) + q\ / e~ * {r2+2ry) / g(x - r\u)du 

Jo Jo 



dr. 



where g(x\y) = if x < 0. The first term accounts for the case where the first seed is V-type, 
and the second term for the case where it is H-type and the process is effectively re-started 
with a different boundary condition having already covered some horizontal distance. Note that 
the ray length probability density function is g{x | 0). 

2. The Moment Generating Function 
We define the moment generating function of the conditional density g(x \ y): 

poo 

M t (y) = / e tx g(x\y)dx, 
Jo 

which, from equation ([T|), satisfies the integral equation: 

roo 

(1) M t (y) = (l-q) (x + y)e-^ x2+2xy Ux 
Jo 

cr+y / poo 



-Ur 2 +2ry) 



e x g{x — r \ u)dx du 



dr. 



To solve this equation we use the following method. We begin by converting equation ([T]) into 
a differential equation whose general solution is a combination of Kummer confluent hypergeo- 
metric functions. Retaining only those parts of the solution which have the correct asymptotic 
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behaviour, it remains to deduce the form of a single arbitrary function of t. This function is 
found by substituting the solution back into the integral equation evaluated at y = 0. 

Commencing with the first part of our method, we note that the second term of ([T|) contains 
the integral: 

/■oo poo 

I e xt g(x - r \ u)dx = / e ( - z+r ^ t g(z \ u)dz 
Jo Jo 

= e rt M t (u). 

We may now simplify the second term using integration by parts: 

rr+y 



e _I( r 2 +2r( j / _ t) 



M t (u)du 



o 



dr 



erfc 



y-t 

V2 



M t (u)du + 



erfc 



u — t 



M t (u)du 



Evaluating the integral in the first term of (pQ), and then combining both terms in together we 
obtain: 



M t (y) = (1 - q) 



1 + y| tei( ,-»)- erfc (^ 

y-t 



erfc 



V2 



M t {u)du + 



erfc 



u — t 



M t (u)du 



This may be reduced to the differential equation: 



dy 2 

which has the general solution: 
(2) M t (y) = 1 + c(t)H q ^ 



d 2 M t (y) . sdMAy) ^ . %r . . 
; -( y -t)—p^-(l-q)M t (y) 



dy 



y-t 

V2 



+ d(t) 1F1 



i-g i (y_ ±_ 

2 '2'^ >/2 



where c(-) and d(-) are arbitrary functions of t. The last term of expression ([2]) is an example 
of a Kummer confluent hypergeometric function. Because all the t derivatives of this function, 
evaluated at t = 0, diverge as y — > oo, we may discard it from the solution on the grounds that 
for all x > 0, liniy^oo g(x \ y) = 0. The second term in the solution, a generalisation of the 
Hermite polynomials, H n (x), to non-integer n, is in fact a combination of Kummer functions 
with the correct asymptotic behaviour [3]. For our purposes, we will need the following integral 
representation of the Hermite function: 

2 v+l ^ r co ^ 2 

e 



(3) 



H v (z) 



TTV\ 

~2J 



Now that we have have the correct general form of the moment generating function, it remains 
to evaluate the constant c(t). To do so, we substitute the solution ([2]) with d(t) = 0, into the 
original integral equation ([T]) evaluated at y = 0. We find that: 



c(t) 



t erfc (^) 



1+2 

-e 2 



In order to find an explicit expression for c{t) we need to evaluate the integral: 



J(t,q) :-- 



erfc 



u — t 



H, 



</-i 



u — t 



du. 
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To do this, we first make the change of variable z = (u — t)/y/2, and then split the integration 
range into two, giving: 



/oo 
erfc(»iT 3 _iO)d. 

= V2 



2" 



£ + §^ - q 2 T{ _ q/2) 



Now consider the integral: 



K(t, 



t 



evic(z)H q -i(z)dz. 



Integration by parts gives: 

K(t,q) = 1 [erfc(*)flg(z)]° * + -4= T e- z2 ^(z)dz. 
2<? QV 71 " V- * 

In order to evaluate the second term, we make use of the integral representation of the Hermite 
function ([3]), finding that: 



rO 99+1 roc fO 

/ e~' HJz)dz = — =■ / e~ u u q / 



2 g+i 
2" +1 



/ 7T?; , 
cos \2zu — ) 

73 



du 



: ^-u 9 1 sin ( | cos ( - (-nq + v^iii) ) du 



\V2 



2"- 



■t cos 



The function K(t,q), in full, is therefore: 



■2q 



r 



+ 



2<? 
2q~^ 



+ 13 ^ 



2 ' 2' 2 



K(t,q) + 



T{-q/2) 21 



4T(1 - (?) q 2 T(-q/2) 



and J(t, g) is: 

J(t,q) = V2 
Putting this all together we have: 

M t (y) = l + 

Wi° *"B*-i (75) 

Setting y = we have an exact expression for the moment generating function for ray lengths 
in the half-rectangular Gilbert model. It is worth noting that our analysis has not lost any 
generality by considering an east moving ray. To obtain the moments of a south moving ray we 



t erfc (^) 
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Moment 


Recurrence 


Exact 


1 


1.81696 


1.81696 


2 


4.64107 


4.64107 


3 


15.57 


15.5701 


4 


65.9719 


65.9721 


5 


342.236 


342.243 



Table 1 . The first five moments of terminal ray length when q = | computed us- 
ing the first 200 coefficients from the Cowan-Maa recurrence relation, compared 
to the exact results from the moment generating function. 



make the switch q f-> 1 — q. Using this expression we find that the first four moments of terminal 
ray length are: 



Ml 



M2 



M3 



I) 



V2T (1 - §) 



r i 



3r(i 



I) 
(*-§) 



+ 2 



v^r (i - 1) 



1 + 2g + 



lj 4 = 8 



l + q + q(l + 2q) 




These have a pleasing form, all being polynomial is the same ratio of Gamma functions. We note 
that the expression [i\ represents a considerable simplification of our previous formula for this 
moment [3]. To verify that our analysis is correct, let us compare these exact expressions to the 
approximate moments calculated using the first 200 coefficients from the Cowan-Ma recurrence 
when q = |. Table [T] shows the results of these calculations, which indicate that our exact 
expressions are correct. 

We conclude by noting that when q = \ the moment generating function has a particularly 
simple form: 



M(t) 



r(4) i*i 



4' 2' 2 



+ V2tr (I) 1 f 1 



4' 2' 2 



3. Concluding comment 

We have determined the exact form of the moment generating function for terminal ray length 
in the rectangular half-Gilbert model. To the authors' knowledge, this model remains the only 
Gilbert-style model for which analytical results exist. 
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